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Abstract 

We consider a convection process in a thin loop. At Ra = Ra'^j. a first transition 
leading to the generation of corner vortices is observed. At higher Ra a coherent large- 
scale flow, which persists for a very long time, sets up. The mean velocity v, mass 
flux m, and the Nusselt number Nu in this flow scale with Ra as v oc m oc Ra^'"^^ and 
Nu oc Ra^'^, respectively. The time evolution of the coherent flow is well described by 
the Landau amplitude equation within a wide range of i?a- variation. The anomalous 
scaling of the mean velocity, found in this work, resembles the one experimentally 
observed in the "hard turbulence" regime of Benard convection. A possible relation 
between the two systems is discussed. 



1 Introduction 



Thermal convection in a Benard cell is one of the classic, well-controlled, systems, on which 
one can test the theoretical understanding of various natural phenomena like fluid insta- 
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bilities, transition, strong turbulence itself and the laws governing heat and mass transfer 
in turbulent flows. Since convection is one of the most commonly occuring phenomena in 
Nature, the ability to describe it is also of great practical interest. The linear stability of 
an infinite fluid layer between two plates heated from below was investigated more than a 
Century ago and the appearence of convection rolls has since been well understood. More 
recent results on transition to turbulence in Benard cells demonstrated a beautiful and in- 
tricate picture of chaos onset from ordered and coherent rolls. Studies of high i?a-number 
turbulence in a convection cell were typically based on the idea that the temperature profile 
averaged over horizontal planes Q{z) = T{x, y, z) differs from an almost-constant- value only 
within thin thermal boundary layers of width 5t- Then, assuming that the upper and lower 
boundary layers do not "communicate" one obtains on dimensional grounds: 

5t ~ Ra~^ (1) 
leading to the scaling of the dimensionless heat flux (Nusselt number, defined below) 

Nu ^ Ra^ (2) 

The accurate experiments on Benard convection, using low temperature helium, con- 
ducted by the Libchaber group, Heslot et al. (1987), Castaing et al. (1989), showed that 
this relation is, in fact, incorrect and instead the Nusselt number at Ra > 10® scales as: 

Nu Ra"" (3) 

with X ~ 0.285, which is close to a; = 2/7. The experiments demostratcd the appearence of 
this scaling as a transition from a random state of the fluid with the heat transfer dominated 
by small-scale vortical motions, observed at 10^ < Ra < 10^. to a new state of turbulence, 
characterized by the onset of a powerful and persistent coherent large scale-flow (boundary 
layer "wind"). This wind, fascilitating strong correlation of the top and bottom boundary 
layers, leads to deviations from the "classic" scaling exponent ^ = 1/3. 

Since the first experiments, Heslot et al. (1987), Castaing et al. (1989), this effect 
was observed in Benard convection with air, helium, water and mercury as working fluids, 
different aspect ratio cells etc, Wu (1991), Wu and Libchaber (1992), Belmontc et al. (1994). 
The experiments revealed that the dependence of the "wind" velocity with Ra is: 

V oc Ra^ (4) 
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with ^ 0.48 — 0.49 differing from the expected free-fall exponent ^ = 0.5. 

The proposed theoretical models mainly dealt with the explanation of the observed 
X — 2/7 exponent, assuming the existence of the large scale flow. Understanding the rea- 
sons for the appearence of coherent flow in strongly turbulent Benard convection was the 
prime motivation of this work. However, the unexpected results, presented below, are of an 
independent interest disregarding their relation to turbulent convection. 

Our motivation was prompted by the following qualitative argument. We consider a 
fluid with Prandtl number, Pr = uq/kq ^ 1, where Uo and k,o stand for viscosity and thermal 
difFusivity. Assume that the flow in the prc-transition ("soft turbulence") state is simply 
a traditional Kolmogorov-like turbulence. We treat the flow as consisting of two parts: a 
viscous boundary layer of width 5 ~ 5^ and the bulk, where the effective transport coefficients 
are estimated as: 

Thus, the bulk can be perceived as a very viscous (large mass) fluid. Setting, for the sake 
of the argument, z/ ^ oo we conclude that the problem of stability of the thin boundary 
layer adjacent to the walls of the convection cell can be decoupled from the "super-stable" 
chaotic flow in the bulk. A somewhat similar situation was considered in the end of the 
sixties by Welander (1967) and Keller (1967). interested in a convection process in a fluid 
contained within a thin loop heated from below. It has been shown that when Ra was large 
enough this system is unstable and a steady clock-wise or counter-clockwise flow sets up. In 
some situations an oscillatory solution was possible. 

The physical explanation of this effect is as follows: The mean bouyancy force in this 
system is approximately, Welander (1967): 

B ^ Agpoa J Tdy (6) 

where A is a cross-sectional area, g is the gravity acceleration, a stands for the thermal 
expansion coefficient, T for a temperature distribution and y is a coordinate in the vertical 
direction. The bouyancy is balanced by a mean friction force F, depending on the ffow rate 
q = vA (where v is the mean velocity) 

F (X q (XV (7) 
when q is not too large. The most important part of the argument is that the buoyancy 
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B = B{q). Model equations, developed in Welander (1967) and Keller (1967), showed that 
B{q) was only weakly g-dependent and, as a consequence, the curves B{q) and F{q) must 
cross at least at one point. That is the reason why a regular flow sets up in this system. 

Thus, the qualitative picture of turbulence in Benard cell, presented above, combined 
with an idea that the low viscosity boundary layer is, in some respect, decoupled from the 
bulk makes the analogy with convection in a thin loop possible. However, if this analogy 
holds, it must explain the experimentally observed anomalous scaling V{Ra), from (1). This 
is a prime goal of this paper. 

We investigated the dynamics of the flow in two-dimensional cells < x,y < L with 
viscosity z/ = z/q < cxd in the interval 6 < \x\,\y\ < L — 6. Outside this interval 1^ = 00 
meaning that v = 0. Experiments showed that the main contribution to the heat transfer 
(more than 50%) came from the wind. That is why at this stage we neglected the heat 
tarnsfer in the inner part of the cell, setting there k, = 0. 

In order to test the sensitivity of the results to the geometry of the system, we have 
performed simulations in three configurations, corresponding to cells of different geometry. 
The first configuration is a cell of L = 1 and S/L = 0.1, the second is identical with the first, 
except for the aspect ratio 6/L = 0.05 instead of 0.1, whereas the third one corresponds 
to a circular geometry with the same diameter L and aspect ratio with the first case; in 
the latter case, the bottom fourth of the circumference of the circle is maintained at a high 
non-dimensional temperature of 1, whereas the top fourth is kept at a low temperature of 0. 
The three convection cells and the boundary conditions used are shown in Fig. |l]l,II,III. 



2 Formulation 



The equations of motion are the Boussinesq equations, outlined below: 



(9v 1 

— + V ■ Vv = -Vp + T + — W 
dt Re 

(8) 

dT 1 
dt RePr 

(9) 

V-v = 

where the non-dimensionalizing velocity scale is U = u/LGr^^"^ and Re = Gr^^"^, where Or is 
the Grasshof number (here, Ra = Gr since Pr = 1). The only non-dimensional parameter in 
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the system (except for shape and aspect ratios) is the Rayleigh number, Ra = PgATL^/ ku. 
The results from simulations we have performed will be reported as function of Ra, and in 
particular as function of r = {Ra — Racr)/Racr, where Ra^r is a critical Rayleigh number 
where transition to large scale mean flow occurs. 

For the time integration of equations (]T0|), we use a fractional step method, in con- 
junction with a mixed explicit /implicit stiffly stable scheme of second order of accuracy in 
time, Karniadakis et al. (1991). A consistent Neumann boundary condition is used for the 
pressure, based on the rotational form of the viscous term, which nearly eliminates split- 
ting errors at solid (Dirichlet) velocity boundaries, Tomboulides et al. (1989). The spatial 
discretization of the resulting Helmholtz equations is performed using two-dimensional Leg- 
endre spectral elements. Patera (1984). The resulting matrices for the numerical solution 
of the two-dimensional Helmholtz equations are solved using preconditioned conjugate gra- 
dient iterative solvers. The resolution can be increased by either increasing the number of 
elements or the order of the interpolants inside each element. In the simulations presented 
here, the resolution was improved mainly by increasing the order of interpolants. Several 
resolution tests, not reported here, were performed and a typical discretization consists of 
40 elements with up to 15 points in each direction per element. In general, because of the 
laminar nature of the flow in the range of Ra numbers investigated, the computational cost 
was not a limiting factor; a typical simulation took few hours on a SGI/RIOOOO workstation. 
The resolution became limiting only in the simulation of very high Ra cases (over 10^ or so) 
described in the last section. 



3 First transition 

At very low Rayleigh numbers, the fluid inside the cell is not in motion. As Ra increases 
over Ra'^^, the flrst transition, corresponding to the appearance of convection rolls, occurs. 
Here, because of the geometry, the convection rolls at the scale / ~ 5 appear only at the 
corners of the domain, and display a double-flip symmetry with respect to the diagonal. 
Isocontours of vorticity, consisting of counter-rotating vortices located at the four corners, 
and for Ra < Racr (corresponding to the second transition), are shown in Figures ^,11, HI. 

This flow conflguration is typical for all Ra'^^. < Ra < Racr- This range of Ra numbers 
was not investigated further, since our interests lie in the study of the large scale mean flow 
which appears for Ra > Racr when the symmetry of the flow disappears. 
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4 Second transition 



As the Ra number increases, a symmetry-breaking, resulting in the appearance of a large 
scale mean flow, occurs. This transition is a linear one and corresponds to a regular bifur- 
cation (or exchange of stability) where the resulting flow is not time-dependent (i.e. the 
crossing eigenvalue has zero frequency). A typical behavior of the total kinetic energy of the 
flow is shown in Fig. ^ which corresponds to case I for Ra = 30, 000 {Racr for this case is 
equal to 28,289). 

Isocontours of the x, y velocity components u, v, temperature, T, and vorticity u, are 
plotted for cases I, II, and III, in Fig. |[ The generation of a large scale mean flow is 
evident when comparing, e.g. the vorticity with Fig. ^. The direction of the flow (clock- or 
counter-clockwise) is random and depends on the initial round-off error disturbances. For 
example, it can be observed from the same Figure, that for case I the large scale flow has a 
clockwise rotation, whereas the opposite is true for the other two cases II, and III. 

The average mass flux (m) and non-dimensionalized heat flux {Nu) for case I are plotted 
in Figures ^and ^, respectively, as functions of r = {Ra — Racr) / Rcicr- The value of Racr for 
case I, was found to be Racr = 28289. As can be observed, the mass flux scales approximately 
as m oc r^''^^ for a large range of Ra numbers (r between 0.1 and 10). However, very close to 
Racr, 1^ oc r°-^ as expected. Landau (1987). The scaling exponent of Nu also approaches 1 
as Ra — i> Racr- In section ^ an analysis is presented in the context of the Landau amplitude 
equation which explains the two types of behavior. It has to be noted that to obtain the 
steady state value of the mass and heat flux for cases close to Racr the equations of motion 
had to be integrated for very long non-dimensional times and steady state results for r below 
0.01 {logijr) = —2) were not performed due to the computational cost involved. This is the 
reason that one can only observe the anomalous 0.45 scaling for case II, since the transition 
in the scaling exponent from 0.5 to 0.45, in this case, occurs below r = 0.01. The variation 
of m and Nu for case II is shown in Figures ^ and ^ respectively; here, Racr = 114648 and 
it can be seen that the r°'^^ scaling is present for the whole range r investigated. The same 
holds for Nu which behaves as r^'^ for 0.01 < r < 1. 

At Ra < Racr, where v = 0, the temperature profile is symmetric with respect to the 
plane a; = and the temperature distribution is a solution to the Laplace equation; this 
solution is very close to a linear profile in y. As flow with an average velocity v = m/5 starts 
to develop, the temperature profile becomes more and more asymmetric and the temperature 
distribution in y— deviates from the almost linear profile at Racr- As expected, as convection 
continues to increase the temperature profile steepens close to y = 1 (for clockwise rotation 
and along the left vertical channel). This steepening is governed by the following equation: 
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with boundary conditions T = at y = and T = 1 at ?/ = 1. The solution to this 
equation, Tg = {exp{vRa^^'^y) — exp{vRa^^'^)^ / (l — exp{vRa^^^)^ , for different values of 
vRa}/^, is compared with temperatures profiles along y (along x = 0.05 which corresponds 
to the middle of the left vertical channel) obtained from the numerical simulations. This 
comparison is shown in Fig. |^ and as can be observed the simple model (p!0D can describe 
the steepening of the temperature layer close to y = 1 from Ra = Racr up to r = 10. 

After this steepening of the temperature profile, the relation rh oc r^'^^ is not valid any 
more. This occurs for values of vRa^/'^ > 50. After this, the mass flux, m, stops increasing 
and its value stays at approximately the same level. On the other hand, a mixing process sets 
up, leading to a homogenization of the temperature along the left and right vertical channels. 
The average temperature in the two channels becomes equal at very high Ra ^ 10®, and the 
long time solution of the equations is not steady any more. It is shown in section |^ that in 
this range of Ra, the thermal energy input can no longer be converted to a coherent large 
scale flow but is almost entirely converted into small scale features and eventually turbulence. 

In order to explore this anomalous scaling m oc oc r°'^^, and to verify that it is not 
due to numerical errors, or due to numerical singularities (at corners), we analyzed case III 
which consists of a perfectly symmetric and non-singular geometry. Here, the value Racr 



was found to be equal to 34, 417. Figures |10| and [TJ show the variation of m, and Nu with 
r, respectively. As can be observed from these figures, the scaling exponents in this case are 
the same as the ones of case I, for both m and Nu. For all cases it was found that m oc r^/^ 
for r < 0.1 or so, and m oc r°'^^ up to r < 10. It was also found that the non-dimensionalized 
heat flux, Nu oc r for approximately r < 0.1, and Nu oc r"'^ for 0.1 < r < 10. 

The results for case III are very close, both qualitatively and quantitatively, to those for 
case I. On the other hand, case II differs significantly from both cases I and III; the rh oc r°'^ 
and Nu oc r range has not been observed. We believe that in this case this range is too 
narrow to be detected numerically. Thus, all we observed was the rh oc r°'^^ and Nu oc r°'^ 
range for the whole range of r investigated. 



5 Scaling of Nu with Ra 

One can use an integral form of the energy equation to obtain an estimate of the scaling 
of Nu number with r. Using the steady state results for all Ra > Racr, we used a control 
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volume which consists of the top half of the convection cell, as shown in Figure 12. The 
energy equation integrated over this control volume is 



/ (n ■ v)Tcis = — ^ / ^ds (11) 
Jcv RePr Jcv on 

The convective heat flux is non-zero only at the lower part of the control volume along 
sides 1 and 2 as shown in Fig. whereas the diffusive heat flux is only non-zero at the 
top boundary T and sides 1 and 2. Also, since the vertical walls at interfaces 1 and 2 are 
insulated, the temperature profiles in the x— direction are very close to constant and therefore 
the x— dependence of T at 1 and 2 can be neglected. One can also add and subtract the 
mean velocity u (which is the same at both 1 and 2 due to incompressibility) and obtain the 
following expression from equation |ll] (by assuming an arbitrary clockwise circulation): 

_,,..,(f,_f,)__L_(|H-f) (12) 
which after non-dimensionalization with the conductive heat flux Qcond becomes: 



u6 (Ti 
1 = — ^ 



Nu-1= \ + I + 1 / I 2^ 1 - 1 



Qcond 




(13) 



The second term, T2, in the right hand side of equation (|13D (inside the square brackets) 
is equal to when r = (conduction regime). When Nu >> 1, this 0{1) term becomes 
negligible. Therefore, when Nu >> 1, one would expect a scaling of Nu — 1 very close to 
the scaling of the first term, Tl, in the right hand side. Our results indicate that v oc r°'^ for 
r < 0.1, and r°'^^ for 0.1 < r < 10. We have also found that the difference (Ti — T2) oc r°'^ up 
to r = 0.1 and then starts leveling off up to r = 10 at which point it starts decreasing, due 



to the mixing process. Figure |13| shows the variation of the logarithm of these two terms, 
as well as the variation of their sum which should be equal to Nu — 1 in a steady flow. It 
can be observed that both terms start increasing as r^'^ at r << 1 and for r > 1 the second 
term becomes negligible. It can also be seen that as long as the flow at long times is steady, 
the sum of these two terms equals Nu — 1 (shown with big open circles) and its scaling is 
also ° for r up to 0.1 and r°'^ after that and up to r = 10. 
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6 Amplitude Equation 



After the second transition at Ra = Racr, a large scale mean flow is generated, the amplitude 
of which increases with Ra. The kinetic energy of this flow, A^, integrated over the whole 
domain was used as an indicator of transition and the Landau amplitude equation was used 
to model this transition. The kinetic energy is governed by the following amplitude 
equation: 

^ = - aA^ (14) 
dt ' ^ ^ 

The solution to this equation is given by the following expression 



A^(t)=7(^^e-(*-') + a) (15) 



where 7 and a are the so-called Landau constants, and Aq = A{t = to)- We found that 
our results can accurately be modelled using equation (^) for a wide range of Ra numbers 
much higher than Racr- It is clear from (ffl) that at steady state A is given by 



A ^ 



and since we know the amplitude A from the DNS, we varied the value of 7 to obtain the 
closest description of our data. Figure n shows the time history of the total kinetic energy 



of the flow, A^, for case III and Ra = 40,000 or r = 0.16. The solid line corresponds to 
numerical simulation results, whereas the dotted line is obtained using equation |l5|; as can 
be observed from the figure the difference is almost negligible. This is true for r up to about 
10 or so, a result which is quite surprising when taking into acount that equation ([T5|) is 
only valid close to transition and more surprisingly that the same equation models data 
demonstrating the anomalous scaling of r°'^^ within this range of r. 

The Landau constants 7 and a were calculated from numerical simulation results and 



are shown in Figures cases I, II, and III, respectively. These figures suggest that for 
cases I and III considered, 7 oc for r — > (as expected for the kinetic energy), whereas a 
is approximately constant in this range and equal to approximately 100; at about r = 0.1 or 
so the scaling exponent for 7 changes to about 0.9 (shown with solid line) and this scaling 
persists up to r between 5 and 10. This is consistent with the results presented in the 
previous section m oc r°'^ for r < 0.1 and m oc r"'^^ after that. The second constant a stays 
almost constant during this transition. 
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For values of r higher than about 5 or so a transition is observed in the scahng of 7 oc r*^-^ 
scahng is no longer valid. This transition also affects the magnitude of the second constant 
a, (which is supposed to be independent of Ra) in the interval of _Ra-variation after the 
second transition, and now this second constant starts to decrease. After this value of r 
equation ([l5| ) no longer describes the data. In fact, for values of r greater than 10 one can 
clearly observe the appearance of other modes which are oscillatory and although damped 
at long times, they do appear in the initial transients. These modes lead to instabilities at 
Ra numbers higher than 10^ for case I, and after that the flow becomes time dependent. 

For case II again the 7 oc dependence was not observed. Instead, the value of 7 was 
found to scale with 7 oc r^'^ for 0.01 < r < 1. Atr>l this relation breaks down, in contrast 
to the rh oc r°'^^ law which is valid up to r = 10 or so. In addition, equation (15) is still 
a good approximation to the time variation of the kinetic energy for at least up to r = 10. 
The second exponent a was found to be very close to 190 for all r up to 1 and to decrease 
monotonically for higher r. 

The most unusual feature of the systems considered above is that the Landau equa- 
tion description, originally proposed for the immediate vicinity of the transition point, is 
remarkably accurate in the entire range < r < 10. 



7 Higher Ra corresponding to r >> 10 

Values of r significantly higher than 10 were investigated only for case I. It was found that 
even for Ra numbers at which the long time solution is steady, oscillatory modes appear in 
the initial transients. In some cases these modes are similar to the pulsating instabilities 
analyzed in Welander (1967), Keller (1967), where the flow may actually reverse its direction 
(clock- or counter-clockwise) before it settles in a steady state. An example is shown in Fig. 
[18| where the time history of the m— velocity component is shown at two different points, one 
in the middle of the lower channel and one in the middle of the upper horizontal channel 
for case I and for Ra = 10^. As can be observed from this figure, the flow reverses once at 
a non-dimensional time of around t = 20 and another time at t = 50 before it approaches a 
steady state after a damped oscillatory behavior. 

As Ra increases to Ra > 10^, the flow becomes unsteady. This unsteadiness is likely the 
result of finite amplitude disturbances. The flow turns abruptly generating the disturbances 
at the corners. However, now the viscosity is not sufficient to damp these disturbances and 
they work as a continuous source of noise which drives the channel flow unsteady. Patera 
and Orszag (1983). Results at these high Ra numbers are shown in terms of vorticity and 
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temperature isocontours in Figure 19 for Ra = 10^ and 10^ respectively. As we can see 
from these Figures, the entire thermal energy input is now transformed into finite amplitude 
vortical structures instead of organized large scale flow. The existence of these structures may 
give rise to three-dimensionality and eventually turbulence via secondary type instabilities, 
Patera and Orszag (1983). 



8 Summary and Conclusions 

We considered a simple system which can mimic the experimentally observed large-scale flow 
generation in Benard convection. This is a two-step process. First, the instability leading to 
a symmetric steady flow pattern happens at Ra = Ra^r- The instability of this pattern leads 
to a symmetry-breaking large-scale flow. Close to Racr the flow rate fn oc r°'^ and Nu oc r 
which is the expected result. 

After the symmetry-breaking the observed w oc rii oc RaP''^^ and Nu oc RaP'^ . The onset 
of this anomalous scaling is correlated with the simultanious modification of the Rayleigh 
number dependence of the coefficients in the Landau amplitude equation, accurately describ- 
ing the data in an extremely wide range of the i?a-number variation. The observed effect 
and numerical values of the exponents seems to be insensitive to the geometric features of 
the system at least for the three cases considered in this paper. At the present time we do 
not fully understand the physical origins of the anomaly. 

The qualitative similarity of the systems considered in this paper with the high Ra num- 
ber Benard convection may explain the experimentally observed large-scale flow generation 
leading to transition from "soft" to "hard" turbulence in a convection cell, Heslot et al. 
(1987), Castaing et al. (1989). This analogy is supported by the similar anomalies in the 
scaling of the "wind" velocity V{Ra) observed in both systems. 
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Figure Captions 



Figure |T]: Geometric configuration of tfie tfiree convection cells investigated 

Figure ^ Formation of convection rolls for the three configurations I) Ra=10,000, II) 
Ra=100,000, III) Ra=30,000. Because this figure was generated by transformation of 
color plots to grey scale, the darkest tone does not correspond to the highest value. 

Figure ^: Total kinetic energy of the flow in time for Ra = 30, 000, case I, a) lin-lin 
plot, b) log-lin plot, which shows the linear theory regime 

Figure ^ a)Isocontours of the velocity in the x— direction u, b) in the direction 
V, c) temperature T and d) vorticity u for cases I, II, and III, and for Ra = 100, 000, 
200, 000 and 200, 000, respectively. The minimum and maximum values of all variables 
are noted in each plot. Because this figure was generated by transformation of color 
plots to grey scale, the darkest tone does not correspond to the highest value. 

Figure ^ Comparison of temperature profiles along y for x = 0.05, for case I with 
equation llOl 



Figure Logarithmic plot of the large scale average mass flux as function of r = 
{Ra — Racr)/ R^cr, for case I 

Figure |^: Logarithmic plot of the Nu number in the large scale flow regime {Ra > 
Racr), as function of r, for case I 

Figure ^: Logarithmic plot of the large scale average mass flux as function of r = 
{Ra — Racr)/ Ro^cr, for case II 

Figure Logarithmic plot of the Nu number in the large scale flow regime {Ra > 
Ra^r), as function of r, for case II 

Figure Logarithmic plot of the large scale average mass flux as function of r = 
{Ra — Racr) / Racri for case III 



Figure |ri|: Logarithmic plot of the Nu number in the large scale flow regime {Ra > 
Racr), as function of r, for case III 



Figure 12: Control volume over which integration is performed 



Figure |T3|: Variation of terms Tl, T2, Nu — 1, and T1+T2 (from equation |T3D with 
r. Solid line is line with slope 0.9 and dotted with 1.0 
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Figure Time history of kinetic energy A^, for case III and r = 0.16 {Ra = 40, 000) 
Solid line is from DNS and dotted line from Landau amplitude equation 



Figure [T5|: Logarithmic plot of the Landau constants 7 and a as functions of r 
{Ra — Rttcr) / Rdcr in the large scale flow regime {Ra > Racr), for case I 

Figure |16|: Logarithmic plot of the Landau constants 7 and a as functions of r 
{Ra — Racr) / Rdcr in the large scale flow regime {Ra > Racr), for case II 

Figure [1^: Logarithmic plot of the Landau constants 7 and a as functions of r 
(-Ra — Racr) / Rdcr in the large scale flow regime {Ra > Racr), for case III 



Figure |T8|: Time history of m— velocity at the points shown, for case I and -Ra = 10 

Figure |T9|: Isocontours of temperature T (a,c) and vorticity u (b,d) for case I and 
-Ra = 10^ and Ra = 10^, respectively. The minimum and maximum values of all 
variables are noted in each plot. Because this figure was generated by transformation 
of color plots to grey scale, the darkest tone does not correspond to the highest value. 
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